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Gamma ray bursts may be blueshifted bundles of the relic radiation 
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A hypothesis is proposed that the gamma-ray bursts (GRBs) may arise by blueshifting the emis¬ 
sion radiation of hydrogen and helium generated during the last scattering epoch. The blueshift 
mechanism is provided by such a Lemaitre - Tolman (L-T) model, in which the bang-time function 
ts(r) is not everywhere constant. Blueshift arises on radial rays that are emitted over regions where 
dts/dr 7 ^ 0. The paper presents an L-T model adapted for this purpose and shows how it accounts 
for the observed properties of the GRBs; some properties are accounted for only qualitatively. 
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I. MOTIVATION AND BACKGROUND 

In the Lemaitre [l| - Tolman Q (L-T) models, in 
which the bang-time function tsir) is not everywhere 
constant, radial light rays emitted close to those points 
of the Big Bang (BB) at which dts/dr ^ 0 display 
blueshifts to later observers (the blueshift is infinite, 
z = —1, on rays emitted exactly at the BB [l-Q)- 

On the other haird, gamma-ray bursts (GRBs) are ob¬ 
served and are believed to originate at large distances 
from our Galaxy, up to several billion light years . The 
question thus arises: could GRBs have been emitted in 
the last scattering epoch, together with the relic radia¬ 
tion now observed as the cosmic microwave background 
(CMB), and then blueshifted to gamma-ray frequencies 
by the mechanism mentioned above? 

For the blueshift mechanism to work, the GRBs would 
have to originate in regions that emerged from a locally 
delayed BB.^ The relic radiation is emitted a finite time 
after the BB, so its blueshift must be bounded from be¬ 
low {z > 2niin > ”!)■ The technical problem to solve is 
this; Is Zmin sufficiently small that, with the free func¬ 
tions of the L-T model suitably chosen, the frequencies 
are blueshifted from the range of the emission spectra of 
hydrogen and helium (the only elements present in large 
amounts during last scattering) to the gamma-ray range 
observed today? The present paper attempts to answer 
this question in the positive - see Sec. m 

Section |TT] provides the most basic information on the 
GRBs. Section Hill is an introduction to the L-T models, 
and Sec. m provides information on light propagation 
in these models. Section |V] presents the current best-fit 
L-T model that reflects the properties of the GRBs, and 
discusses improvements in the model needed to achieve a 
full quantitative fit. Conclusions are summarized in Sec. 
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^ Regions where the BB occurred earlier than in the background 
generate shell crossing singularities iQ]. 


ED The Appendix presents the data needed to replicate 
the numerical calculations. 


II. BASIC FACTS ABOUT THE GRBS 

The following properties of the GRBs need to be ac¬ 
counted for 

(1) Their frequencies extend from ~ 0.24 x 10^® 
Hz to i^^max ~ 1-25 X 10^^ Hz (see also Ref. my 

(2) They last from less than a second to a few minutes. 

(3) They are probably focussed into narrow jets. 

(4) A GRB is sometimes followed by a longer-lived and 
fainter “afterglow” at larger wavelengths. 

(5) Nearly all GRBs come from very large distances, 
from over 10® to several billion light years. 

Currently, there is no generally accepted explanation 
of origins of the GRBs. There exist only attempts at 
explanation by known astrophysical phenomena such as 
gravitational collapse to a black hole, a supernova explo¬ 
sion or a collision of ultra-dense neutron stars 

The model presented in Sec. Elaccounts quantitatively 
for the lower limit in property (1) and for (5), quali¬ 
tatively for (3-4), and is not in contradiction with (2). 
References to these properties will be marked there by 
bullets ® . To achieve a quantitative agreement, more 
elaborate fitting will be needed. 


III. THE L-T MODELS 

The metric of the L-T models is: 

ds^ = df - , dr^ r)(di?2 -b sin^ ddip^), 

1 -1- 2E{r) 

(3.1) 


^ Converted from keV to Hz by: leV = 1.6 X 10 = h X 0.24 X 

10^® Hz, where h = 6.626 X J s 
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where E{r) is an arbitrary function. The source in the 
Einstein equations is dust; its (geodesic) velocity field is 

= (5“o. (3.2) 


Because of the assumption p = 0 built into this model, it 
is inadequate before the last-scattering epoch. 

The function R{t,r) is determined by 

R/ = 2E{r) + 2M{r)/R, (3.3) 


M(r) being another arbitrary function; we neglect the 
cosmological constant. We will consider only the models 
with > 0 and E > Q. The solution of (13.31) is then; 


R{t,r) 


sinh rj — rj 


M 


(cosh 77 — 1 ), 


(2B)3/2 

M 




(3.4) 


where ts {r) is one more arbitrary function; the BB occurs 
at t = tsir). The mass density is 


up = 


2M^r 
R^R^r ’ 


def SttG 
K = - 


The r-coordinate is chosen so 


that @ 


M = Mqt^, 


(3.5) 


(3.6) 


and Mq = 1 (kept in formulae for dimensional clarity). 

The units used in numerical calculations were intro¬ 
duced and justified in Ref. [l^- Taking Q 

1 pc = 3.086 X 10^3 km, 1 y = 3.156 x 10^ s, (3.7) 


the numerical length unit (NLU) and the numerical time 
unit (NTU) are defined as follows: 


rescale the affine parameter A so that 

k\to) = -1. (4.3) 

Then, with being given by (ixa . we have 

l + z = -k^. (4.4) 

For nonradial rays, on which C ^ 0, the last term in (14.11) 
will go to infinity when i? —0. Thus, at the BB 

lim |fc*| = 00 . (4.5) 

Equations (HU) and (14.41) imply that 2 ; —>■ 00 at the BB 
on all nonradial rays, in agreement with Ref. [^. 

For radial rays, z—>■—Iasi?—>-0 when dts/dr ^ 0 
at the intersection of the ray with the BB, and 2 —>• 00 
when dts/dr = 0 The result 2 ; = — 1 implies 

infinite blueshift for all observers. Rays emitted close to, 
but not right at the BB will acquire a finite blueshift 
that can be overcompensated by later-acquired redshifts 
if the observation is carried out sufficiently far to the 
future from the emission point. 

In consequence of (14.51) . blueshifts may arise only on 
radial rays. Thus, for an observer re-directing her tele¬ 
scope away from the direction to the center of the ra¬ 
diation source, the transition from blueshift to redshift 
would occur abruptly - if the real Universe were exactly 
modelled by the L-T geometry. In reality, the changeover 
from blueshift to redshift can be expected to occur in a 
finite, but short time. This would account for the short¬ 
livedness of the GRBs (property (2) in Sec. HIl), their 
narrow-jet appearance (property (3)) and for their “af¬ 
terglows” (property (4)) - see Sec. |V]for details. 

In the following, past-directed radial rays are dealt 
with, on which C = 0. Using 6ID, the equations to 
be integrated numerically are: 


I NTU = I NLU = 9.8 x I0^° y. (3.8) 


IV. LIGHT RAYS IN AN L-T MODEL 


The geodesic null vector fields in (13.ip obey [l^ 


{kf 


i?,.^ (fc”)' 
1 -f 2E 



with C being constant along a geodesic. 
The general formula for redshift is 0 




(4.1) 


(4.2) 


where k^ is an affinely parametrised vector field tangent 
to a light ray connecting the source and the observer, 
both comoving with the cosmic medium with the four- 
velocity u“. The subscript “e” means “at the emission 
event”, “o” means “at the observation event”. We will 
consider past-directed rays, on which k* < 0. We can 


dE 

’ dA “ ’ 


(4.6) 

^/1 + 2E 

dr 

(4.7) 

'' ]J 

R,r “ ’ 

_ — pr 

dA ’ 


with the initial condition (14.31) . The sign in (14.71) is -)- on 
past-inward rays and — on past-outward rays. 

In a general L-T model with E ^ 0 we have Q 

R,r = (^-^)i? + 'k(t,r)i?,*, (4.8) 

As can be seen from (HU) and dUl), R,tr = 0 is the locus 
of extrema of redshift along radial rays; we call it the 
maximum-redshift hypersurface (MRH). 

Consider a radial ray proceeding to the past from an 
initial point that lies later than the MRH. The redshift 
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z on it increases from 0 to a maximum, achieved at the 
MRH. Further down the ray, 2 decreases. If the ray could 
continue to the BB, z would decrease to —1. However, the 
L-T model ceases to apply at the last-scattering hyper¬ 
surface (LSH). Can z become sufficiently negative before 
the ray crosses the LSH for shifting the optical frequen¬ 
cies to the gamma-ray range? It is shown in Sec. IVlthat 
this is indeed possible when the functions Eir) and tsir) 
are suitably chosen, and the observer is put in the right 
place at the right time. 


V. FITTING THE L-T MODEL TO 
OBSERVATIONS AND MEASUREMENTS 

The mass density at the instant of last scattering in 
the standard ACDM model is [l^ 

Kpis 88 X 10^ (NLU)"2. (5.1) 

We assume that the recombination occurs at the same 
density also in an inhomogeneous Universe. The density 
along a ray is calculated using (j3.5|) , and the value of 2 ; at 
the moment when p = p\s emerges from (14.411 and (14.611 . 

The GRBs cannot arise by blueshifting the whole 
black-body spectrum of the relic radiation to the gamma- 
ray range. First, the spectra of the GRBs do not have 
the black-body forms (example: Ref. 0 )- Second, the 
intensity of a GRB created in this way would exceed the 
observed ones by tens of orders of magnitude.^ So, if the 
GRBs arise by blueshifting the relic radiation, then the 
different frequencies have to be blueshifted individually. 

To shift the lowest frequency of the hydrogen emission 
spectrum 0 , J^Hmin = 4.054 X 10 ^^ Hz (corresponding 
to the wavelength of 7400 nm) to the lowest observed fre- 
miency of the gamma-ray bursts, y^rain ~ 0.24 x 10 ^® Hz 
0, the blueshift 1 + 1.667 x 10 “® is needed. To shift 

the frequency of the most intense line in the hydrogen 
spectrum, 2.1876 x 10 ^® Hz (corresponding to 656.2852 
nm) to the same minimum gamma-ray frequency, 1 + z 
1.9 X 10 “^ is needed. To shift the maximum measured 
hydrogen emission frequency, i^Hmax = 3.2 x 10 ^® Hz (cor¬ 
responding to 93.782 nm) to the maximum recorded cos¬ 
mic gamma-ray frequency, ~ 1-25 x 10 ^^ Hz Q, 

I + z fv 2.56 X 10 “® is needed.^ 

The configuration shown in Fig. [T] was obtained by 
fitting the functions E{r) and tsir) and the position of 
the observer by trial and error so as to make 1 -|- z as 
close to zero as possible. The current best result is 

1 + = 1.23007568 x 10“® (5.2) 


^ The Planck formula and the current CMB spectrum Ea imply 
1.295 x 10® W/(cm^x sr x Hz) for the maximum intensity of 
the black-body spectrum blueshifted from LSH by 1 -|- ^ = 10“®. 
Observed GRBs have intensities below 10“^“^ W/(cm® Hz) [l3l . 

^ The most intense helium emission lines have wavelengths be¬ 
tween 388 nm and 846 nm, i.e. within the range of the hydrogen 
spectrum 1201 . 


between the LSH and now. This accounts for the lower 
end of frequencies in ® Property (1). This is the mini¬ 
mum z within the chosen class of BB profiles; other pro¬ 
file classes may possibly lead to smaller 2 . 



FIG. 1: The blueshifted ray, the bang-time profile and the 
profile of the maximum-redshift hypersurface in the L-T 
model described in the text. See the Appendix for the pa¬ 
rameters of ts(r). The MRH is symmetric around the r — 0 
line, but its left part is suppressed in this graph. 

For simplicity, E was chosen the same as in a Fried¬ 
mann model, with 

2E/r^ - k = -0.4. (5.3) 

A general E would have —k replaced by (—fe -|- E{r)), 
where J'(O) = 0, but otherwise is arbitrary Q. This 
would provide further parameters for fine-tuning. 

The BB profile consists of a spherically symmetric 
hump surrounding the center of symmetry out to a finite 
distance; further away from the center ts is constant, 
and so the geometry is Friedmannian. The present time 
is t = 0, and the flat part of ts was chosen at 

t =^tBf = -0.13945554689046649 NTU 

ft! -13.67 X 10® years; (5.4) 

this is the asymptotic value of ts in the L-T model that 
imitates accelerated expansion using nonconstant ts 0 . 

A radial cross-section through the hump in the BB is 
shown in Fig. [21 It consists of two ellipse arcs (see Fig. 
[3] in the Appendix) that are connected by a straight line 
segment in the neigbourhood of the point, where the full 
ellipses would be tangent to each other. The lower ellipse 
arc is tangent to the flat part of the ts profile. This shape 
is determined by five parameters: four semi-axes of the 
two ellipses, and the tilt of the straight segment (see the 
Appendix for their values). 

The inset in Fig. |2] shows the neighbourhood of the 
maximum of ts(r). The ray passes Ate = 0.125 x 10“^ 
NTU = 1.225 X 10^ years over this maximum. The Ate 
is also an adjustable parameter. The blueshift at the 
observer position is sensitive to the value of Ate in certain 
ranges. For example. Ate = 0.126395 x 10“® NTU results 
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FIG. 2: Main panel: Closeup view of the neighbourhood of 
the BB in Fig. [T] Inset: Closeup view of the neighbourhood 
of the BB maximum - see text for a comment. 


in 1 + z ft! 3.9 X 10-5, while Ate = 0.1264 x lO'^ NTU 
results in 1 + z ft 4.5 x 10“^. Thus, a change in Ate by 
5 X 10“® NTU ft 49 years causes that 1 + z goes up by the 
factor 0.87 x 10^ - enough for the observed radiation to 
drift from the gamma-range into the ultraviolet [2l|. The 
49 years is not a sufficiently short time to account for the 
abrupt end of a GRB quantitatively,® but the qualitative 
effect of afterglow is here, and thereby ® Property (4) 
is qualitatively accounted for. 

But the abrupt beginning and end of a GRB and the 
afterglow should rather be associated with the center of 
the radiation source going into and off the line of sight 
than with the passage of time at the observer (see Sec. 
CyI). The model predicts a discontinuous jump from red- 
shift to blueshift and back to redshift for an arbitrarily 
small change in the direction of observation (which qual¬ 
itatively accounts for ® Property (3)). In reality such 
a change might be provided, for example, by the orbital 
motion of the Earth. Then, the blueshifted ray would 
stay within the observer’s field of view only briefly, and 
this would implicitly account for ® Property (2). 

If the center of the BB hump would stay in the ob¬ 
server’s line of sight all the time, then she would see 
the gamma radiation persisting for nearly the whole 
Ate = 1.225 X 10^ years. (But this is a property only of 
the concrete BB profile of Fig. [2] This does not exclude 
the existence of profiles providing shorter viewing times.) 
It would show up abruptly, because the ray emitted near 
the top of the hump would have large redshift. 

The model of Fig. [1] cannot explain the highest- 
frequency end of the GRB spectrum Q by means of 
blueshifting from the hydrogen emission range. This does 
not yet imply that the high-frequency GRBs must arise 
by a different mechanism: it is still an open question 


® Moreover, the accuracy of numerical calculations for this paper 
was insufficient to capture time intervals of the order of hours. 


whether a finer-tuned L-T model could do the job. 

A cosmological model that would account for the mul¬ 
titude of observed GRBs should be imagined as a Fried¬ 
mann background containing many humps like the one in 
Fig. [2J of different shapes, different spatial extents and 
different heights above the flat part of tsir), placed at 
different comoving positions. 

Redshift is used in astronomy as the measure of dis¬ 
tance. But blueshifting renders redshift - distance rela¬ 
tions multivalued [l^. No operational method of deter¬ 
mining the distance is known when blueshifts are present. 
The distance from the observer to the source of the ray 
in Fig. [T] can be estimated in two ways: 

1. The intersection of the ray with the LSH occurs 
ft 13.3764 X 10® years ago, so by conventional accounting 
the source would be 13.3764 x 10® light years from the 
observer. This is 6.83 x 10® years later than the BB of 
the ACDM model, whose present age is 

r= 13.819 X 10® y = 0.141 NTU. (5.5) 

2. The sum of the values of the r-coordinate of the ob¬ 
server (ro = 0.41946) and of the intersection of the ray 
with the LSH (rpsHi = 0.09841) is = 0.51787. The 
re determines the active gravitational mass contained 
within the r = Ve sphere centered at the observer by 
M = MoVe^. In the ACDM model, the same rg (cor¬ 
responding to the same mass) is reached by the present 
observer’s past light cone at t = te = —0.1279 NTU, 
i.e. ft 1.25 X 10® years ago,® which implies the distance 
ft 1.25 X 10® light years. This number is within the range 
of distances inferred from observations of the GRBs Q, 
and thereby ® Property (5) is accounted for. 

The most obvious possibility to improve the model is 
to manipulate the shape and size of the BB hump seen 
in Fig. [U for example, by putting more parameters into 
it. This may result in the rays staying before the MRH 
for a longer segment of the affine parameter, and thus in 
smaller observed values of 1 -I- z. 

The second possibility is to model the BB hump by the 
Szekeres metric Mm 0- It contains the L-T model 
as a subcase, but in general has no symmetry. Therefore, 
it must be verified what directions in the Szekeres class 
of metrics go over into the L-T radial directions, and 
whether blueshifts appear on them. A Szekeres model 
might produce a different time-profile of the observed 
frequency, and a better quantitative agreement with the 
observed properties of the GRBs. 


VI. CONCLUSIONS 

The model proposed here deals satisfactorily with the 
low end of the GRB frequency range (property (1) in Sec. 


This calculation was based on the numerical results of Ref. 1231 . 
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El) and with the low end of the distance range (property 
(5)). It also deals qualitatively with poperties (3) and 
(4). Re (3): an arbitrarily small misalignment between 
the line of sight and the direction to the center of the 
radiation source causes the gamma-ray impulse to disap¬ 
pear. Re (4): the afterglow is present in the model of 
Sec. El but its duration does not agree with observa¬ 
tions. The only property not explicitly accounted for is 
(2). In order to account for it, one should calculate the 
change of observed frequency induced by a small change 
in the direction of observation. This is possible to do in 
the present model, but requires numerical calculations of 
a much higher precision. 

Using a BB profile with more parameters, and mod¬ 
elling the BB hump using a Szekeres rather than the L-T 
metric may lead to a full quantitative agreement between 
the model and the observations. 


Appendix A: The BB profile 

The hump in the BB profile consists of two ellipse arcs 
and of a straight line segment joining them that passes 
through the point where the full ellipses would touch each 
other - see Fig. The values of the parameters in Fig. 


Hare Bq = 0.01, Bi = 0.09, Aq = 0.016, Ai = 0.018, 
Xq = 0.000859. The other ones are uniquely determined 
by these. The units for Aq and Bq are NTU; the other 
parameters are dimensionless. The values of Aq and Bq 
imply for the time difference between the maximum of 
tB and its flat part 0.026 NTU = 2.548 x 10® years (0.18 
of the age of the Universe given by (15.51) 1. 



FIG. 3: Parameters of the bang-time profile (drawn not to 
scale for better readability). See text for the actual values. 
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